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ABSTRACT 

We present radio coiRiumim mapping, optical imaging and .spectroscopy of the newly discovered 
double-peaked broad-lined AGN WISE J233237,05— 505(543.5 at redshift z — 0,3447. Tins source ex- 
hibits an FR.-I and FIt-11 hybrkl-inorpliolog> r t characterized by bright core, jet* and Doppler-boosted 
lobe structures in ATCA continuum maps at L5, 5.6 t and 0 GHz. Unlike most FR-II objects, 
W2332— 5050 is basted by a disk-like galaxy. The core has a projected 5" linear radio feature that 
is perpendicular to the curved primary jet, hint ing at unusual and complex activity within the inner 
25 kpe. The multi-epoch optical-near-IR photometric measurements indicate significant variability 
over a 3 21) year baseline from the AGN component. Gemini- South optical data shows an unusual 
double-peaked emission-line features: the centroids of the broad-lined components of Ho and Hd are 
bhicshifted with respect to the narrow lines and host galaxy by ~ 3800 kms” 1 * We examine possible 
cases which involve single or double supermuss ive black hoicks in the system, and discuss required 
future investigations to disentangle the mystery nature of this system. 

Subject handings: Galaxies: individual: WISE .1233237.05—505043.5; Galaxies: jets; Galaxies: nuclei: 
Galaxies: active; Galaxies; interactions; Radio continuum: galaxies 


I. INTRODUCTION 

The discovery of supermassive black holes (SMBHs) 
with masses greater than 10 billion solar masses in two 
nearby galaxies (McConnell et at. 2011, 2012) highlights 
the need to understand how SMBHs grow and influence 
galaxy evolution. Two commonly accepted growth pro- 
cesses arc act ret ion from the interstellar medium (ISM) 
mid SMBH mergers (Alexander le fflckOK 21)12 L The lat- 
ter scenario is thought to play a more dominant role for 
SMBH growth in the early universe, though binary AGN 
seen at late stages in the merging of bulge- dominated 
galaxies should still be relatively common in the nearby 
Universe, However, despite significant efforts searching 
for resolved binary AGN (e.g_, Komnssa et ah 2003; Ro- 
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driguez et ah 2t)(Mi; Civano et al. 201(1; Gram et nl, 2010; 
Fit et ah 2011a, 2012; ComtTford cl ah 2012; Barrows et 
ah 2012) and late-stage (unresolved) merging SMBHs (.lu 
et al. 2013: Shuu et al* 2013), such systems appear to be 
very rare. Two potential signatures of merging SMBHs 
art 1 from their complicated radio morphology (see Banin- 
garte k: Shapiro 201 1 ) and unusual spectroscopic features 
such as velocity offsets between broad and narrow spec- 
tral line components (see review by Popovic 2012; Ko- 
nmssu 2012). 

^BMH binaries arc the inevitable result of galaxy for- 
mation. In a major galaxy merger event, the SMBHs 
hosted in two large merging galaxies dynamically inter- 
act as the parent galaxies coalesce and form a final bulge. 
SMBH merging can be considered in three stages (see 
Merritt k Milosavlje vie 2005; Uogdimovicel ah 2008, and 
references therein). In the iirst stage, the SMBHs dissi- 
pate their angular momentum t [trough dynamical fric- 
lion with surrounding stars, thus reducing their orbits* 
Tliis process is efficient near the nuclear region of t he host 
galaxy where the stellar density is high, bur stalls at radii 
of ~ 0.01 If) pc due to the depletion of stars. SMBHs 
spend tin* majority of their lifetime at this dose sepa- 
ration (Begelinaii et al, HisO: Yu 20112). However, if the 
system contains a gaseous accret ion disk, the binary orbit 
could decay further via interaction between the SMBH 
binary and the gas (Escula et ah 2004, 2005; Dotti et ah 
2U07 S 200!)). This is the so-called “intermediate stage” , 
lasting for a few to a few dozen Myr* In the final stage, 
when two SMBHs are brought close enough, the binary 
orbital angular momentum dissipates efficiently via grav- 
itational radiation and the SMBHs proceed to final coa- 
lescence. 

Multi-jet radio morphology is potential strong evidence 
for SMBH binaries or SMBH mergers. In tin- ease of 
3C 75 (z 0.0232), four radio jets are launched from 
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a pair of arm-ting AGNs separated by S kpe (Owen et 
al 1985; Yokosawn it I none 1985), Another case, the 
radio galaxy 3C GOB (z = 0.0215), shows zig-zag pat- 
terns of compact radio jets associated with an SMBH 
binary, separated by 4 x 10 _;| pc in a 1.1 year orbital 
period (Sudan ct id, 2003). Igiiehi et al. (2010) calcu- 
lated the decay time of 3C G6B orbit is about 500 yr due 
to gravitational radiation, although a pulsar timing ex- 
periment did not detect the expected gravitational wave 
from tills binary to a 95% confidence level (Jeuet et al. 
2004), From their short lifetimes, SMBH binary systems 
with small orbital separations and short orbital periods 
are expec ted to be very rare. X -shaped radio galaxies 
(XRGs; Leahy k Parma 1002), a subclass of Fanaroif- 
Riley Class II (FR-II; Fnuaroff k Riley 1074) objects, 
contain two misaligned pairs of Centro-symmetric radio 
lobes originating from an elliptical galaxy host . Usually 
one pair of more extended lobes with low surface bright- 
ness is oriented at a large angle with respect to a pair of 
less extended, high surface brightness lobes. Their radio 
morphology has been suggested as a sign from rapidly 
changing SMBH spin, possibly tracing an SMBH merg- 
ing event (Dennett -Thorpe et til. 2002; Merritt k Ekers 
2002; Gopal- Krishna el al. 2003). Other proposed expla- 
nations for XRGs include back-flowing material from jets 
deflected off ISM, outflow from over-pressured cocoons, 
ur interactions of jets with stellar and mtergalactic shells 
(Leahy *1* Williams 1984; Worrall et al. 1995; Capetti 
et al 2002; Hodgrs-Kluck et al 2010). These models 
have 1>< tested in X-ray and optical spectroscopic mor- 
phological studies (e.g* SuripaUi k Subrahmanyan 2009; 
Cheung et al. 2009: Loudt et al 2010; Hodges- K luck et 
al 2010)* A rec ent review of observational evidence ar- 
gues that none of the proposed models can explain all of 
the observed properties of XRGs (Gopai-Krlslma et al 
2010 ). 

In optical spectroscopy, host galaxies with broad emis- 
sion lines significantly Doppler-shifted from the corre- 
sponding narrow lines have been interpreted as the con- 
sequence of a late stage SMBH merger, or alternatively, 
due to a kick received by the broad- lined SMBH result- 
ing from anisotropic gravitational emission (Gonzalez et 
al 2007; Civano ct al. 2010), Theoretical models pre- 
dict a maximum “kick" velocity of ~ 4000 kins -1 for a 
merger inspiraling on a qiiasi-circulur orbit (CaiupaiielH 
*4 al 2007}* In certain unusual orbit, configurations such 
as hyperbolic encounters, the velocity displacement AT 
between narrow lines (NLs) and broad lines (BLs) can 
reach 10,000 km s _I (Healyet al 2009). Observatlonally, 
however, large BL-NL velocity offsets are not commonly 
seen* In the past few years, studies that mine the Sloan 
Digital Sky Survey (SDSS) spectroscopic archive have 
identified only 16 objects whose broad Baliner lines are 
offset from the narrow lines by > 3500 km s" 1 (Komussa 
ct al. 2003; Boroson k Lauer 2009: Shields et al 2009; 
Sleitihardt et al. 2012: Tsulmnutza el al 2011; Bracleous 
et al. 2012), In addition, two alternate models have also 
been proposed to explain this unusual spectral configura- 
tion in these objects: (1) a Uno-of-sighl superposition of 
two unrelated AGN (c.g., Shields et al 2009), and (2) an 
asymmetric Kepleriau accretion disk (Chen k Halpern 
1989). 

The accretion disk model, which requires a centrally- 


illuminated, geometrically thin, and optically thick ac- 
cretion disk, has successfully explained the double- 
peaked BL emitters (“double-peaked emitters” or DPEs; 
Bomiing cl it!. 2007; Shields et al 2009; Eracleoiis et al 
2009), It is expected that accretion disks should natu- 
rally become circular after a period of time (Reis 1988; 
Caimizzo et nl 1990 ), and for this reason the circular, 
Kepleriau disk model of Chen k Halpern (1989) lias been 
used with success on many DPEs. However, an elliptical 
accretion disk (Eracleons et al, 1995) might arise if an in- 
stability is introduced in the disk, such as tidal disruption 
of a star or an orbiting body (such as a second black hole) 
which might allow an elliptical accretion disk to persist 
for some time (Syor Clarke 1992). The disk profiles 
are typically centered around B aimer emission line po- 
sitions presumably because of photoionization and elec- 
tron scattering in the high-density medium of the disk. 
There Is typically a peak on the red side and a peak on 
t he blue side of the narrow emission lines. The specific 
shape of the line profile* is dependent on geometrical fac- 
tors such as the disk inclination, inner and outer radii, 
line broadening, and surface emissivity. The blue peaks 
can sometimes be stronger t han the red peaks because of 
modest relativistic beaming. 

Other configurations of an accretion disk can also pro- 
duce similar double-peaked BL features. For example, 
asymmetric BL profiles can raise from non-imifonnly dis- 
tributed giis such as that in long-lived eccentric accre- 
tion disks (Ogilvie 2091)* The warped accretion disk, 
produced by radiation-driven Instability (Maloney et at* 
199G), can generate the double-peaked line emission pro- 
file (\Vu et al 2008), In addition, a one-arm spiral excited 
by self-gravity in the accretion disk can also create asym- 
metry in the line profiles (Storehi-Bergtnaim et a l 2003). 
These configurations, although all capable of producing 
double-peaked BL features, predict different detailed line 
profiles and variability, and can be tested by long-term 
spectroscopic monitoring (Gezari et al 21)07: Lewis et 
at. 2010). A few other complex scenarios involving two 
distinct orbiting SMBHs have also been proposed to ex- 
plain a few individual DPEs with unusual profiles (Tang 
k Grind lay 2009; Barrows et ah 2011)* 

In tins paper, we present the newly discov- 
ered double- peaked emitter WISE *1233237.05—505643,5 
( W2332— 5056 hereafter) at s = 0.3447* W2332-5056 
has a radio luminosity of 2.0 x 10 2r> W Hz“* str -1 at 4.85 
GHz, comparable toother FR-II galaxies (FanaroflTl: RK 
ley 1974, assuming a spectral index a — —0.8 between 
4.85 GHz and 178 MHz). The radio continuum morphol- 
ogy of this object is complex, exhibiting both FR-I and 
FR-II characteristics in addition ro a winding primary jet 
with multiple structures near the core, including a linear 
structure perpendicular to the primary jet. The optical 
spectrum of W2332— 5056 exhibits broad emission lint's 
(Ha and H r d) blneshifted by 3800 km s^ 1 with respect 
to the corresponding narrow lines. The combination of 
rhese features makes W2332— 5056 a prime candidate for 
an SMBH merger system. 

Wo describe the multi-band observations of 
W 2332 -5050 at mid-infrared (rnid-IR), radio, op- 
tical, and X-ray wavelengths in §2. The detailed 
results comprise §3* In [j3.7, %ve interpret our results 
and discuss various interpretations of this unusual 
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object. Wc summarize our work in §4, Throughout 
this paper, we assume a concordance cosmology with 
Hq = TOkmiT 1 Mpc~\ tt m = 0.3, and fl A = 0,7. For 
this cosmology, 1" at z = 0.3447 subtends 4.9 kpc. 

2. SOURCE SELECTION AND OBSERVATIONS 

W2332— 5056 was identified based on its unusual mid- 
IR colors from the NASA Wide- field Infrared Survey Ex- 
plorer ( WISE; Wright et til. 2010). WISE mapped the 
entire sky in four mid-IR bands (3.4, 4.6, 12 and 22 pm, 
or Wl, W2, W3, and W4, respectively), achieving typical 
point-source sensitivities of 0.08, 0.3 1, 1 and G mJy (5a) 
in the four bands, respectively. In the short-wavelength 
bands, W2332-5056 has a relatively red |3.4]-[4.6j color 
of 0.9 mag, similar to normal AGN (Stern ct at. 2012) 
and occupies the same color-color space as other SMBH 
merger/rccoiling candidates discussed in §1. However, 
the [4.6]- [12] color is 2.0 mag T significantly bluer than 
the typical AGN seen by WISE (Jarrctt ct al, 2011). It 
is also a strong radio emitter, having been previously de- 
tected at 4.85 GHz in the PMN Survey (Gregory et. aL 
1094; Wright ct al. 1994) and at 843 MHz in the SUMSS 
survey (Matidi ct at. 2003). We obtained higher spatial 
resolution radio follow-up observations using the Aus- 
tralian Telescope Compact Array (ATCA), and optica) 
imaging and spectroscopy with the Magellan, SOAR and 
Gemini South telescopes. Details of the observations are 
given below. 

W2332-5056 is 8.7 ajcinmutes SE of the background 
z ~ 0.5707 galaxy cluster SPT-CL J2331-5051, which 
was discovered by the South Pole Telescope (SPT) 
Sunyacv-ZePdovich survey (Vandcrliudc ct al. 2010). As 
a result, it has been observed with Magellan optical 
broadband imaging (Inamori Magellan Areal Camera 
and Spectrograph, or 1MACS; Dressier et al. 201 1), the 
Spitzer I RAC camera, and the Chandra ACIS-I camera 
as part of SPT cluster follow-up programs. These obser- 
vations arc detailed below, and aperture photometry is 
presented in Tiible \ . 

2.1. Radio Observations 

The radio continuum observations were carried out at 
ATCA using the Compact Array Broadband Backend 
(CABB) with 2x2 GHz bandwidth. W2332-5056 was 
mapped at 1.5 GHz in the GC configuration on UT 2010 
July 1G-17. The 5.6 GHz and 9.0 GHz data were obtained 
simultaneously on UT 2010 July 18-10 in the GC config- 
uration (maximum baseline ^ 6.0 km), and on UT 2010 
August 10 in the hybrid H1G8 configuration (with 5 an- 
tennas closely positioned within 200 m of each other, and 
the 6th antenna at 4.5 km distance). Observations were 
taken in snap-shot mode over I2h (6C) and 6h (HIGH) 
periods. A total of 25 minutes integration was obtained 
at 1.5 GHz, and 166 minutes at 9.0 GHz and 5.G GHz. 
PKS 1934-63 was used as the flux calibrator while PKS 
2311-452 (for 5.6 GHz and 9.0 GHz) and PKS 2333-528 
(all three bands) were used for phase calibration. The 
largest angular scales (#las) visible under these config- 
urations arc 5ff> at 1.5 GHz, 3f9 at 5,6 GHz, and 213 at 
9.0 GHz. 

The data reduction was done using miriad. The final 
robust- weigh ted images are shown in Figure 1 (a)-(c). 
The primary b emu at the highest frequency, 9.0 GHz, 


is > G\ significantly larger than the 1" field of interest 
for W2332-5056, and hence we did not apply a primary 
beam correction. The rms of the 1.5 GHz map is 0.62 
m.Ty beam -1 with a beam size of I3?12 x 4"80, PA = 
6*. The (u,v) data at 5.6 GHz and 9.0 GHz from the 
two configurations were combined, resulting in an rnis 
noise* of 29 pJy beam^ 1 and 37 j*Jy beam^ 1 in the final 
synthesized beams of 2?85 x T/51, PA = 13° and 1780 
x 0787, PA = 8°, respectively. The results of the inte- 
grated flux measurements using aifs are fisted in Tabic 
t Because the 0 ca$ for each band is larger than the an- 
gular extent of W2332“'5056, we anticipate negligible or 
minima] missing flux in the ATCA interferometric mea- 
surements. 

2,2. Mid-IR Observations 

Tkble 1 presents mid-IR photometry of W2332-5056 
from the WISE All-Sky Release Atlas Catalog (Cutri et 
al. 2012). The observations of this field were carried out 
on UT 2010 May 19-20. A total of 12 single frames in 
each WISE band are coadded. The spatial resolution of 
WISE is -6" in Wl, W2, and W3, and -12" in W4. 
W2332— 5056 is well detected in Wl— W3, with signal- 
to-noise ratio (SNR) > 15, and marginally detected (4 
<r) in W4. The mid-IR fluxes are listed in Thble 1 . 

In addition, W2332-5056 was observed in Wl and 
W2 during the WISE post-cryogenic survey on UT 2010 
November 13, 16, 17, and 18. A total of 19 frames were 
acquired, spread across 3 days due to the satellite ma- 
neuvering to avoid the Moon during its Sun-synchronous 
polar orbit. These 2nd epoch measurements allow us to 
test the variability of the AGN in W2332-5056 at mid- 
IR wavelengths, as discussed in §3,-1. 

Finally, due its projected proximity to SPT-CL 
J 2331 -5051, W2332— 5056 was also observed in [3.6] 
during Spitzer Cyclc-9 (ID: 60099, PI: Mark Brodwin) 
on UT 2009 November 26. The pipeline-processed data 
(PBCD; shown in panel (d) of Figure 1), retrieved from 
the Spitzer Archive, provides another mid-IR epoch for 
comparison. 

2.3, Optical and Near-IR Observations 
2.3. L Optical and Near-IR Imaging 

Optical photometry of W2332-5056 was obtained with 
Magellan*! M ACS in the SDSS ^/iV-bands as part of 
SPT cluster follow-up efforts (Vandcrlinde et al. 2010). 
The observations were taken on UT 2008 November 3- 
4 with median seeing of ^ 078, and the 5ht magnitude 
limits in AB system are ~ 24.8, 24.8, 24.4, and 23.4 for 
</, i* % and respectively. The details of the Magellan 
observations and data reduction are described in High et 
ah (2010). 

To explore variability, we obtained archival photomet- 
ric data of W2332— 5056 from several major sky surveys. 
The optical photometry from SuperCGSMOS Sky Sur- 
vey (Humbly et ul. 2001) in UKST Bj- r /?-, and /-bands 
has a photometric uncertainty of — 0.3 mag. Near-IR 
imaging from the Two Micron All-Sky Survey (2MASS; 
SkruLskie et ul. 2006) was obtained on UT 1999 Octo- 
ber 18 mid provided 5 ~ 9 a detections in the J-, H-, 
and AV bands. Detailed information of these sky survey 
observations is presented in Table I . 

In order to study the morphology of W2332-505G, we 
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TABLE L 

MULTJ-WaVK LENGTH PHOTOMETRY OK W2332-5056 


Band R.A. (J2O00.0) Dec). (J2000.0) Wavelength Instrument Flux Density Epoch (UT) Notes 


Hard X-ray 

23 h 32 m 37* * 

Soft X-ray 

*** 

FUV 

23 b 32 m 36?89 

NUV 


UKST-B, 

23 h 32 B “3r06 

UKST-R 


UKST-/ 


9' 

r* 

i' 

z' 

J 

23 h 32 m 3ff96 

23 h 32 m 3r;06 

H 


Ks 

*** 

I RAC Channel 1 

23 b 32 m 37?12 

W1 

23 b 32 m 37705 

W2 


\V3 


W4 

AKARI-Sftm 
AKARI -I8»m 
IRAS- 60/jm 
IRAS-lOOum 
150 GHz 


9.0 GHz 

23 h 32 m 3?;06 

5.6 GHz 

23 h 32 m 3r05 

1.5 GHz 

23 b 32 m 3(£86 

4.65 GHz 

23 h 32™ 37^05 

843 MHz 

23 h 32 m 3ff62 


-50° 56'43" 

2-10 keV 


0.5-2 keV 

-50° 56' 42'/ 8 

1516 A 


2267 A 

-50°56'43'/2 

4400 A 


6600 A 

-** 

8000 A 

-50 c 56'43'/44 

4770 A 

-** 

6231 A 

*** 

7625 A 


9134 A 

-50°56'43'/3 

1.235 |im 

*** 

1.662 pm 

«** 

2.159 pm 

-50° 56' 44'/ 1 

3.6 pm 

-50° 56' 43!' 5 

3.4 jum 

... 

4.6 pm 

... 

12 pm 


22 /j m 

9 pm 
IS pm 
60 
100 
2 mm 

-50° 56' 4 3'/ 6 

3.3 cm 

-50 D 56'43'/3 

5.4 cm 

-50° 56'41'/3 

20 cm 

-50°S6'43'/5 

6. IS cm 

-50° 56' 42" 6 

35.6 cm 


Chnndra/AClS 

53^2 njy“ 
63l}® t>Jy a 

*.* 

GALBX 

6±3 piJy 
7±2 /iJy 

Su perCOS MOS a 

26±6 fijy 
90 ±20 niy 
80±20 /jJy 

Magellan/ IMACS b 

85.1±0 5 *iJy 
182±1 ^Jy 
254±1 fijy 
432±1 >iJy 

2MASS 

0.37 ±0.06 mJy 
0.42±0,09 m 3y 

*.. 

0.79±0.09 mJy 

SpUrer/IRAC 

1.3U±0.002 mJy 

WISE 

1.31±0.03 mJy 
1.70±0,04 mJy 
2.0±0.1 mJy 
4 ±1.0 mJy 

AKARI 

< 0.03 J/ 7 

< 0.05 Jy* 

IRAS 

< 0.5 Jy 

< 1.5 Jy 

SPT Survey 4 

< 2.6 mJy f 

ATCA 

37.2±0.5 mJy" 
48.6±0.3 mJy 1 ' 
103±2 mjy r 

PMN Survey^ 

47 ±9 mJy 

SUMS5* 

175±5 mJy 


2000 August 12* 30 

1 2003 Aug, 10 
J - 2005 Sep. 06 
1975 July 16 
1992 November 14 
1991 August 03 
2003 November 3*4 


1999 October IS 


2009 November 26 

2010 May 19-20 


— 

fwhm = 5 ;: 

- 

FWHM = 5" 

* 

FWHM = if; 

- 

FWHM = 3'( 

- 

FWHM = 1' 

2010 July 18-39 

*i.as = 2.'3 


6 las » 3 '9 

2010 July 16-17 

«las - 5f5 


FWHM = 45 

— 

FWHM = 4/1 


* Comspmidintf to the observed /{0*5-2k*V) = 2,321 £m x 10“ ''i'recm -3 . and /{ 2 - lOkeV) = 1*01+2:*? x 10" 17 wfiem -2 s‘* . 
h Humbly at ni. (2001) 

c Umtt. 

d Joaquin Vieira, private communication. 
c Integrated flux within 0f5 rftdiu.v 

* Gri’R^q' cl nL (JfM)4); Wrifilit ct ul. (1994) 

* Maudi ct *1. (2UO*i) 


also acquired high quality riband imaging {A — 6 165 A} 
witli ff/65 seeing using the 4.1m Southern Astrophysics! 
Research (SOAR) Telescope with the SOAR Optical Im- 
ager (SOI; Walker et ah 2003) on UT 2011 November 22* 
The basic image processing and reductions were carried 
out with I RAF* This r'-band image is shown in panel 
(e) of Figure l* 

2.3*2. Optical Spectroscopy 

The optical spectroscopic observations of W2332— 5956 
were carried out with the Gemini Multi-Object Spectro- 
graph (CMOS; Hook et ah 2004) on the Gemini South 
telescope (Gemini-S) on UT 2011 May 11 and UT 2011 
November 28* We obtained two 10-minute integrations 
observations with the I*'5 slit and the B6G0/520.0 dis- 
perser in May 201 l t covering 3800A to G7G0A with re- 
solving power R = A/AA - 1700* At z = 0*3447, the H/3 
line is near the red end of the wavelength coverage. In 
November 2011, in order to cover H0 and Ha together 
and confirm the broad line feature observed ear- 
lier, we revisited W2332-5056 with Gcmini-S/GMOS, 
with the 1"5 slit and the R400/705J) disperser* The 
wavelength coverage ranged from 6200 A to 0100 A with 


R — 1900. The spectra were processed using IRAF* 
The white dwarf GJ 318 was used for flux calibration in 
May 2011, and the DA white dwarf LTT32I8 was used in 
November 2011* The final combined spectrum is shown 
in Figure 2* 

2*4* X-Ray Observations 

W2332-5056 was serendipitously observed twice with 
the ACIS-I camera in programs to map the galaxy clus- 
ter SPT-CL J233 1-5051* The total exposure time was 
35*0 ks, split between a guaranteed time observation of 
29*0 ks on UT 2009 August 12 (ObsID 9333; FI Garmire) 
and a guest observation of 0*0 ks on UT 2009 August 30 
(ObsID 11738; PI Mohr)* The X-ray data for the cluster 
were reported in Andersson et ul* (2011). The data were 
analysed following standard procedures using the Chan- 
dra Interactive Analysis of Observations (Cl AO; V4.2) 
software. We identified “good time intervals” for the 
observations, yielding a total effective exposure time of 
30*0 ks for W2332— 5G56. 

3. RESULTS AND ANALYSES 
3.1. Radio Continuum Morphology 
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Fin. L \V2332~5Q56 imaging at radio, near- JR, and optical wavelengths. Images are in (a) 1.5 GHz continuum (beam size: 13/1 x -1/3, 
PA = 6°); (6) 5.6 GHz continuum (2/85 x P/51, PA - 13°) ; (c) 9,0 GHz continuum (l/SO x 0/87, PA = 8*); (d) Spitzer IRA C 3.6 fim 
(FWHM ** 1/7), and (e) r'-baod (seeing = 0/65) In the radio maps* (a*c), the colors represent the flux density from peak in each map 
to 3tf level, and the beam sizes of maps are illustrated with blue ellipses in the lower-right corner of each panel. The contour levels are at 
2 H * f| x 5 o, where n is non- negative Integers, and la are 0 47 mJy/beam at 1,5 GHz, 0.026 mJy/beam at 5.6 GHz, and 0,03*1 mJy/beam at 
9,0 GHz. In i lie 3.6 /mi (d) and r'-band images (er), the colors show the flux density peak from background levels on a logarithmic scale. 
The location of the 9 GHz radio core is marked with a ■*+*' in fa), (d), and (e). The final panel (f) shows a composite mid-IR and radio 
view of W2332 — 5056 (blue: Spitzer JR AC 3,6 pm; green: 9,0 GHz; red: 5.6 GHz) in the matched resolution {2/9 x 1/7, PA s 15°). 


The ATCA maps at 1.5 GHz (20 cm) t 5.0 GHz (5.4 
cm), and 9.0 GHz (3.3 cm) arc shown in panels (a)-(c) of 
Figure i. These maps clearly show a mixed morphology 
between FR-1 and FR-I1 structures that comprise the 
W2332— 5056 system. In the lowest resolution map, at 
1.5 GHz (3"x 4'/8), three distinct peaks arc observed, 
including the core and two unevenly bright spots ai both 
ends of the renter. The central peak is offs el by 3/0 from 
the location of the U VST -detected host galaxy (referred 
to as the core), which indicates significant unresolved 
20-eni emission from the extended structure north of the 
core. At shorter wavelengths, the jet is well resolved and 
exhibits a complex morphology, described next. 

In the 9.0 GHz map (1"80 x [f/87, PA = 8°)* we iden- 
tify six distinguishable components, marked as A F in 
Figure 3: A) a core (coincident with the host galaxy), a 
pair of “primary” jets, including B) a south-east (SE) jet 
with PA of 120° and C) the counter distorted north-west 
(N\Y) jet, D) a linear structure at a nearly orthogonal 
PA of 20\ E) a luminous hot spot, and F) an arc ( lose to 
the core. The fluxes of these radio structures are listed 
in Table 2. 

At the position of the mid-lR and optical peak, a radio 
core is apparent. The core has a rt datively flat spectrum: 
a spectral index of a — —0.18 ± 0.07 (F^ a */*) from 5.6 
GHz to 9.0 GHz in the maps of matched beam of 2^9 
x 1"6 at PA = 15°, The deconvolved size of the core is 
< if ! 3 at 9.0 GHz, based on single Gaussian lit ting of the 
uniformly weighted map. 

Originating from the core, the primary jet extends 
about 14" to the north-west , punctuated by a typical FR- 
II galaxy Doppler- boosted hot spot that is much brighter 


than the core. In contrast to the linear jet 0'B") to the 
south-east, the north-west jet shows an asymmetric zig- 
zag morphology. The deviation is Hr from the axis of 
the south-east jet, while the deviation of 30 is shown 
on the opposite side. Moreover, there appear to be dis- 
tinct differences across the radio bands, notably closer to 
the core, better revealed in Figure 1(f). The jet is more 
curved at 5,6 GHz on the south-east side, and relatively 
st might er on the north-west, spatially offset from the 9.0 
GHz emission. Tin 1 nature of tills offset is under investi- 
gation with new, deeper radio observations that we are 
currently pursuing. 

The faint linear structure “D”, pointing to the center 
of the core, extends 5" 7 n from the central core, is vis- 
ible at 5.0 GHz, and notably at 9.0 GHz (appearing as 
•green 1 in Figure' 1(f)), extending over > 2 beams. The 
emission at 1.5 GHz is also extended in the same direr- 
tion indicated by the elongation of the core in the low 
resolution 1.5 GHz map. The co-alignment is shown in 
both 9.0 GHz and 5.0 GHz of the radio core and the 
resolved host in optica! imaging (white/magenta in the 
figure). This structure has signakto- noise > 30. It is at 
a PA of 20°, which is distinct from the PA of the elon- 
gated ATCA beams at 5.6 GHz and 9,0 GHz (13° and 
8°, respectively). Moreover, it does not mated) the un- 
cleaned dirty beam structure of a relatively weak central 
core (peak- to- peak, ratio between core and component 
"D" is ^ 8). Structure "F” has an unusual shape resem- 
bling a curved, arc- like structure, possibly the result of 
a 3-dimension radio structure seen in projection. 

The primary jets extend to 70 kpc on both sides at 
a PA of I20\ To our detection limit, there is no clear 
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FlC. 2. Combined GMOS-S optical spectrum of W2332-5056 from 3800 A to 9100A. The narrow lines are marked at red* hi ft : — 0,344/. 
The sign *©* labels the atmosphere absorption features. The scaled broad Ha and H3 lines arc highlighted in the inset panel, indicating 
velocities relative to the narrow line redshift (shown by the dotted vertical line). The dashed lines around 4715A, 5675A, 7155A, and 8130 A 
represent the CCD gaps of CMOS'S. The signs U V' show the unidentified narrow- band features. 


TABLE 2 

Radio Continuum Flux Densities of Resolved Radio Components 


Component 

ID 

^1 

(mJv) 

Fs.gciu 

(mJy) 

F<i iicili 
(mJy) 

0 5,a 

n y.t> 

Core 

(A) 

1411 

5.0101 

4 610.1 

-0. 1810.07 

Primary Jet (SE) 

(B) 

34 ±1 

11.610.2 

9.510.3 

-o.4i±aos 

Primary Jet (NW) 

(C) 

5711 

27.810.2 

18.110,3 

-0.9010,03 

Linear Structure 

(D) 

— 

1.510.1 

1.110.1 

-0.7010.15 

Hot Spot 

(E) 

*— 

21.210.1 

14,010.2 

-0.8710.03 

Arc- like Structure 

(F) 

— 

— 

1.8310.08* 

— 


Note. — Component IDs (sn nml crilunm) tu*- marked in Figure \. fV ,j ClJl 
mid ft> uciii nr*' pwasurrd in On* I I ruji[h i 279 x i7fl, PA — 15 3 ). 
Tin 1 a, as definM in ft. x y'\ is tin- spectral index between 9.0 GHz and 5.6 

CHfc. 

Measured in in Up of bruin — L'Hfl X 07S7, PA ~ S c . 


jut-like structure that connects directly io the core mi 
the axis of the primary jet in the central 12 kpc (2'/5). 
The western jet exhibits complex behavior which may 
indicate interact ions wit li the interstellar medium (ISM) 
and mtergalacfie medium (IGM), The structure of the 
eastern counter jet is more linear, with a maximum PA 
change of < 10 e . The pitch angle along the jet readies 
20 a in the 5.6 G&l nvUvivus it is within 5° at 

<L0 GHz. If the driving mechanism that produces this 
"zigzagging morphology in the jets Is periodic, the es- 
timated period is on t he order of 0.25 Myr based cm the 


length between off-axis peaks of the serpentine primary 
jet. 

The asymmetric brightness of the jets suggests that 
the western side of the primary jet is on the side near 
to us. This is because jets can be brighter (fainter) by 
a factor of D :i if they are moving toward (away from) 
the observer (discussed in Meier et ah 2001)* where D = 
[F ( 1 t |eosfl)]“* is the Doppler factor for a relativistic 
jet at Lorentz factor F For given velocity v, mid the angle 
between t In* axis of the jet and t he observer's line of sight 
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FIG. 3. 9.0 Gib continuum map with the key radio features 

illustrated: (A) core, (B) primary south-east (SE) jet T (C) primary 
counter north-west (NW) jet, (D) linear structure, (E) hot spot, 
and (F) arc-like structure. 

0. However, 0 for the primary jets is likely to be > 8G Q 
{e.g- the maximum for jets hi the face-on plane) based 
on the relatively small brightness difference between the 
primary jets on cadi side, 

3.2. Spectroscopy: Broad Ha and H3 Emission Lines 

The integrated fluxes of identified lines tire listed in 
the Table 3. Tike presence of the high ionization line 
[NeV]342fiA, as weU as the ratio of |Q 111)5007 A to nar- 
row line H(3 indicate the AGN nature of W2332— 5056. 
Additional evidence for the AGN is provided by broad 
Ho and H 3 line widths, as shown hi Figure 2. Both 
have asymmetric profiles with respect to the NL. 
The full width at half maximum (FWHM) of Ha is 
~ 11,000 kins -1 , and the full width at quarter maximum 
(FYVQM) is ^ 16,000 kms -1 . For Hd, the peak of the 
BL component is centered at — 6453A, or 3800 kins -1 
blueshifted with respect to the systematic host galaxy 
redshift of 0,3447, consistently measured by both the nar- 
row emission lines as well as the CaHK absorption, Simi- 
larly, the peak of the BL Ha is at 8730 A, or 3200 kms -1 
blueshifted with respect to the host galaxy systematic 
redshift. Objects with asymmetric, very broad Balnier 
lines whose peaks arc? shifted from the corresponding nar- 
row lines are usually referred to as double-peaked emit- 
ters (DPEs). We note a few narrow emission features 
(marked us in tin* inset of Figure 2) straddling Hn t 
unlikely to be [N II] or [S II}, which could potentially be 
associated with Ha components shifted by -1700 and 
+4000 kins” 1 with respect to the galaxy systematic ve- 
locity* The different blueshifts for the broad Ha and Hd 
are the result of contamination from these mi-identified 
components. FYoin the difference between the scaled Ha 
and Hd profiles, we detect no significant broad Fe II emis- 
sion near Hd wit liln ± 10,000 kms -1 . We do not. set 1 clear 
evidence for very broad line regions in other spectral fea- 
t tires. 


f-T 



o 


3t,0 36.D 33:32:37,0 36.0 30oO 

Right sBceealon 


3,2. 1 , Pu i u •} n e t art zaii a tt wit h Git ass i a n Profit as 

We first attempted to parameterize the BL components 
with two Gaussian profiles without, excluding narrow 


TABLE 3 

Optical Link Flcxksof \Y 2332- 5056 


Line 

Redbhift 

Flux 

erg cm 2 s’ 1 ) 


Nev] A 3.126 

0.3439 

7.0±0.6 


0 ri[ A 37 27 

03443 

13.4±0.4 


Neill] A3369 

0.3440 

6,3 ±0.4 


O ill] A-1363 

— 

1.2±0.3 n 


Hell A-1686 

— . 

0.S±0.3“ 

H 3 



NL 

0.3447 

2.12±0.2 

BL-blue 


G4,l±0.8 1, 

BL-re<i 

— 

71.8±1.0* 

JO Jill A4959 

0.3449 

U.2±0,2 

pm]A5007 

0.3449 

■12.2±0.2 

O l] AG3G0 

0.3458 

4.4±0.7 

Ho 



NL 

0.3447 

12.3±0,6 

BL- blue 1 

— 

■136±3 il 

BL-red 

— * 

265±2 i; 

[N ll] A6583 

0.3448 

10.3 ±0.6 


11 Mnrglnnl rjoioctimj. 

Inti'graM flux from tho v.vro velocity to tbe blue 
mul of tin* BL compouuut, 

*' luitgrrttt'ij flux from tin- awro velocity to the re<l 
mil of tin* BL comjMUH'ne 


emission lines in order to measure the FHWM and fluxes 
of all lines. The two BL Gaussian profiles and one NL 
component are used for both Ha and Hd in the model, 
allowing to fit independently on Hd mid Ha* The narrow 
absorption feature's (between rest-frame 4866.4 4908* 2 A 
and 6529.3 6552,4 A) are excluded from the fitting pro- 
cess. For the (O m]495f)A and 5007 A lines, one BL Gaus- 
sian component and one NL component are assumed, 
while the |N 1 1)6583 Ais fit. to a single NL Gaussian com- 
ponent, The redshift and llnevvidth of the [0 Ill]4959A 
and [0 IH]5007A lines an* assumed to be equal, with a 
fixed flux ratio F [om ^ Qi]7 / F\oui] 4mj = 3. The results 
(Xm-i — L85, lit >f = 3607) are shown in Figure 1 and 
Table 4. 

The blue BL components are centered at < —4300 
kms -1 with respect to the narrow line, while the red 
component peaks at > +53011 km.r 1 . Both blue and 
ml BL components have FWHM of > 8800 kms -1 . The 
flux of t he blue components is ^ 1.5 2.3 times that of the 
red components. The fitting residuals shown in the lower 
panels of Figure l indicate that the model represents the 
data near both Ha and Hd fairly well, except the blocked 
region and the very red wing of the Hd (AV > 11,000 
kms -1 ). 

3.2.2, Spectroscopic Modeling - the Casa of a Thin , 
Elliptical Accretion Disk 

in the literature, the DPE feature, such as we see in 
W2332— 5056, has often been interpreted with a thin ac- 
cretion disk surrounding a single AGN {Banning el ah 
2007; Shields vi aL 2000; E radix ms ct ah 2000). Wc 1 
test the disk-emitter scenario by fitting the Ha and Hd 
profiles with models based on line emission from an 
elliptical accretion disk (Eraeleous et aL 1095). The 
asymmetric blueshifted broad line features in both Ha 
and Hd in this model are due to Doppler boosted line 
emission from the disk. All NL features (Ha, Hd, 
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Fig. 4. Observed H3 aud Ha broad line* ovcrplotted with Hie best-fit Cau»ian models of BLs plus XLs, The two Gausstau BL profiles 
in the model are fit independently to both of H3 aud Ha The redshift and Unewidth of |0 UIJ4959A and [O II 1)5007 A are assumed to be 
equal, with a fixed fiux ratio F|oiii]soor/^(OtitH(i5i» = 3. The shaded regions which cover the narrow absorption features and wavelengths 
outside the spectral window are excluded in the model fitting. The blue and red dashed lines represent the blueshifted and redsbifted BL 
components of Ha and H3 f respectively. The magenta lines indicate the BL components of (O MI] 4959 A and [O HI]5G07A, The dark green 
dotted lines show NL components of H3 t [O II 1)4959 A. [O 111)5007 A, Ha, and [N IIJ65S4A lines. The green solid line is the sum of all 
components. The lower two panels show the residuals from the fitting. 


TABLE 4 

Emission Link Pa it am etk its 


Line 

AV 

(km s“ 1 ) 

FWHM 
(km=“ 1 ) 

Flux 

(10” 16 erg cm- r ) 

Ha 

NL 

0±4Q 

390±50 

9±1 

BL-blutr 

—44 20 ±50 

918G±30 

519±2 

EL-rcd 

±5330 ±130 

91 10±6G 

229 ±3 

H3 

XL 

— 70±30 

36G±6Q 

1.6±0.2 

BL-blue 

-4270±9Q 

S790±10 

74±1 

BL-red 

+SS60±130 

S600±20 

4S±1 

[Ollll A5007 

XL 

+70±ia 

480±10 

33±1 

BL 

+50±50 

1300±120 

8-1 

IN ii] A 6583 

NL 

— 10±30 

34G+6Q 

6-1 


Note. The results of CatUfsI&n fitting of tin* emission line* iwr 
shown here. Th<* second coin nm, AV, represent the velocity olf- 
Met with rmpfict to the rest frame wavelength of th* 1 line at tec Is hi ft 
of 0,3447, For OAch emission line, one BL component and tme NL 
component b included, exerpt Iht and Hd ILr which two BL compit- 
nents are employed, The fO lii]495f>A (tint shown in the table) and 
[D 111)5007 A lines are forced to have the smile mid FWHM, with 
a 1:3 line finx ratio. 

(0 HI] 4959 A, [O mJfifJOTA, [N Il]6584A} and an unex- 
pected ‘ absor pt ion-like ' region redward of NL H 3 (6535 
9600 A ) are excluded from the model fitting. The eccen- 
tricity, e, tingle of perhistron, o t inclination angle, t, and 
inner /outer radii, r- in / r out , are free parameters for fir- 
t ing. Eac h model lit Ha and H i wit h a radius-dependent 
power-law emissivity e(r) - tu x r"‘ { for the disk, where r 
is the radius of the emission region from the center, mid 
q is the power law index. Two q indices are used ;is in- 
dependent parameters, f/p fl and f/u^, for the emissivities 


of Ha and H3> respectively. Theses emissivity assump- 
tions for the photoionination models (Collin-Souffrm A' 
Dumont 1989; Dumout k Gulihi-Souffrin l!J90a f b) were 
developed by Chen A' Halpern ( 19x9), and are extensively 
used by others, including Enickoits fc Halpern (2003). 

The modeling was performed by fir ting a grid of free 
parameters and adopting the set of parameters that 
yields the lowest reduced chi-squared value. We tested 
elliptical disk models for eccentricities of e — 0.0 0.8 

{show'll in Figure along with I he best lit circular model). 
The high eccentricity e = 0*6 model provides a statisti- 
cally significant improvement to explain the BL profiles 
(\n-ii = 0.88, dof = 1887), although models with e > (1,2 
generally provide good fit [\t*d - 1-02). The best-fit pa- 
rameters of this model are as follows: e = (Ui, t/hc = 2,4, 
qwa = 3,2, i — 12% FWHM local broadening a — 3300 
kins' 1 , rd"St r = 240 r G , and iffi* = 4500 r G , where r G 
is the Schwarxschild radius of the central black hole. 

This model wit li a thin disk and offset BL regions 
can reasonably explain the Ho and H5 line profiles in 
W2332— 5056* other than an apparent "absorption” fea- 
ture at AV ^ 1000 kms“ 3 redsbifted with respect to H/T 
This absorption- like profile is not observed for IIq, likely 
because it overlaps with [Nh] 0583A emission at the sys- 
tematic redshift of the galaxy. The nature of this feature 
is yet unknown. 

3.3, A - ray Emission 

The AGN hosted by W2332-505G is strongly detected 
in Chandra ACI3-I data, with a total of 1790 source 
counts in the energy range 0,3 - TkeV. It. is the bright- 
est X-ray source in the field. In the deeper observation 
(ObsID 9333), our target is 9f3 off-axis, where the Chan - 
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FlG. 5. Observed H.3 and Ho broad lines overplotted with elliptical dUk models* for five eccentricities with all other parameters the 
same: e = 0.0 (blue), c = 0.2 (cyan), e = 0.4 (red), e = 0 6 (green) and e = 0.6 (dark green). Models are fit simultaneously on both Hd 
and Ha with an assumption of a radius-dependeiit power-law emissivity. See text for details. The shaded regions are excluded in the model 
fitting. The e = 0.6 model yields the best fit, though the broad line profiles can be fit relatively well by the e — 0.2, 0.4, and O.S models. 


dm/ ACIS 50% encircled energy radius is 4" at 1.1 kcV, 
We extracted the source distorted by the off-axis PSF us- 
ing a relatively large 12" radius circle. Background was 
from a 30" circular source free region. Response matrices 
and effective areas were then determined for the target 
position, and the spectrum was binned to have at least 
20 counts per spectral bin. The resulting spectrum is 
analyzed using XSPEC (VI 2.0,0). 

W 23 32— 5056 appears spatially extended in tin- Chan- 
dra data, but the distorted off-axis PSF complicates the 
interpretntiou of this spatial structure. We used the 
Cl AO tiusk arestoke to perform PSF deconvolution on 
the deeper image (ObsID 0333) with a PSF simulated by 
CHART In the energy range of 0.3 7 KeV. Two weak, ex- 
tended X-ray regions ^ 47 5 from the central X-ray peak 
at PA of ~ 30 G and ^ 200° are marginally detected: to- 
gether they contains ^ 0.0% of the X-ray photon counts* 

Using the standard definition of hardness ratio, HR = 
(H - 3) /{II + S)> where II is the hard counts between 
2 - UJkeV and S is the soft counts between 0*5— 2keV, 
the derived HR = —0.05 is typical of exiragalaetie X- 
ray sources (e.g., Stern et al. 2002). Our best absorbed 
power-law fit (\ 2 d = 1.00/70dof) to the spectrum has 
an observed photon index P = 1.51 i 0,08, a total ab- 
sorbing column of A'h = 5*31 x 10 21 cm” 2 , an ob- 
served soft X-ray flux of /(0*5 - 2keV) = 2,32tgjg x 
urgan^s” 1 , and an observed hard X-ray finx of 
/( 2 - IDkeV) — l.OliJjY x 10" 12 erg enr 2 s“ l . The 
corresponding intrinsic, or unobscured X-ray luminosi- 
ties are L( 0*5 - 2keV) - L72+SSc x 1044 ergs" 1 and 
L{ 2- lOkcV) — x IQ 44 ergs" 1 . 

3,4. Variability 

W2332— 5056 was observed in Spitzer [3,6] in 20011 and 
in WISE [3.4] between 2010 2011, The median of the 
[3.4] single-frame fluxes is 1,30±0,05 mJy while the [3.0] 


flux from a single measurement is l,311±Q.t)02 mJy; the 
anticipated color difference between IRAC [3.0] - WISE 
[3,4] is less than 0.1 magnitude 1 \ No significant mid- PR 
variation is seen. 

At optical wavelengths, the redshifted H/? spectra were 
obtained twice in 2011, separated by 6 mouths. There 
is no obvious variation in H 3 between these two epochs* 
On longer time spans (>3 years) the central AGN does 
show significant fluctuations at optical wavelengths. The 
residual in the different nil r'-band imaging (Figure 0) in- 
dicates a ~ 20%, total Ilux decrease from 200 k November 
to 2011 November, 

3,5. Optical Imaging Morphology 

W2332—5()5fl is clearly resolved in the SOAR riband 
image (see Figure le and Figure 7: (7/65 seeing)* We use 
GALFIT 3*0 (Peng et ah 2002, 2010) to examine its mor- 
phology* Inspired by evidence of at least one BL AGN 
in the system, we use a combination of Seisie profiles 
and PSFs to model the host galaxy and unresolved AGN 
components, respectively. The results are shown in the 
right panels of Figure 7. The single Sersic model, sin- 
gle PSF model, and dual PSF models can not explain 
the observed morphology well: all of them leave signifi- 
cant residual Ilux (> 3%) within 478 radius of the core. 
Modeling the image with a Sersic profile plus a single 
fret 1 PSF, the PSF preferentially selects an off-center op- 
tical bright spot, but. leaves significant, central residual 
flux (^ 7%). The best fit model ^ dof = 

40388) requires a single Sersic profile accompanied with 
two PSF profiles which match the observed brightness 
distribution to < 0,5%. In this model, the Sersic com- 
ponent, primary PSF component, and secondary PSF 
components contribute 33.7%, 61.8%, mid 4.5% of the 
total flux, respectively, 

* 1 bribed on Figure 4 of 1st tp : / , w . se2 - 1 j ii : : cs . :h . * du / do as/ 
reloiioa/a! l8ky/oipsup/aoc:6_3Q .html 
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LTT 201 1 November 22 
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Fit*. G* — Optica) t* images of W2332— 5056 from two different epochs (left two panels) after convolving to a matched 1"2 FWHM PSF, 
and the residua) difference image (right pane))* The shaded region illustrates the region in the CCD gap in the 2011 observation. This 
region is excluded in the analysis All 5 a sources, including the extended structure of W2332— 5056, are successfully subtracted down to 
the < level. The residual component of W2332-5056, coincident with the nucleus, has a FWHM of l'/2, indicating that it is unresolved. 
The peak flux of the residual is ^ 140 o of the difference image, and contains about 18% of the flux of the peak in November 2008 which 
is brighter than that in 20 Ih 
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FlG* 7.— GALFIT morphological models and residuals images of W2332— 5056 in r' band. The PSF and the original image (taken on 
UT 2012 November 22 with seeing of 0^65) are shown on the left. The other eight panels, in square- root scale, show models and residual 
images (14" on a side). Note that the free single PSF component in the ^rsic ■+■ 1 PSF” model (the second column from the right) selects 
an off-center optical bright spot instead of the central peak. 


The Sersic profile (of the host galaxy) in the best- fitted 
model has an axis ratio of 0*64 at a PA of 4&°* It con- 
tributes ~ 50 /iJy (or AB = 19,64) with ail effective 
radius of l'/3G (6*7 kpc), and has a Sersic index of 0.92, 
close to an exponential disk index n = I* This com- 
ponent has — -21*6 (AB mag) before ^-correction, 
slightly brighter than the Milky Way. The luminosity 
is 2*5 x lO 1 ^ Lq, similar to our SED model of the hast 
galaxy (blue line in Figure H). The primary PSF com- 
ponent is centered 0! 1 6 offset (^ HOD pc in projection) 
from the Sersic center* This displacement, 25% of the 
PSF FWHM, appears significant, but could be the result 
of an imperfect PSF generated from unresolved sources 
in the field. 

The secondary unresolved component is 1*0 (4.9 kpc 
in projection) offset from the Sersic center. Tills compo- 
nent does not show obvious variability* It is too bright 
(A/^ = —19*4) to be a star cluster inside the host galaxy 
or a satellite dwarf elliptical galaxy, though it could be 
a second merging, compact galaxy. It is not likely to be 


stellar object unless it is a foreground star along the line 
of sight* The position angle to the primary compact ob- 
ject is 49°, or ^ 30° with respect to the secondary radio 
jet, making it unlikely to be a bright spot associated with 
an optical jet of the central AGN. Further understanding 
of this component will require sensitive, higher resolution 
imaging* 

The brightest visual companion at slightly greater di- 
sauce (r > 4", ^ 20 kpc) has nv is 21, C (AB mag), 
which corresponds to ^ 6 x 1(P L© in r f band after fe- 
correction assuming the redshift of W2332— 5056* There 
are no obvious disturbed morphological features brighter 
than 1/5 of the luminosity of W2332-5056 to suggest 
that the host galaxy is an carly/iutermediate-stage ma- 
jor merger .system* 

3*6. SED Modeling 

We model the broad-band SED of W2332— 5056 from 
the UV to the mkl-IR with a combination of host galaxy 
and AGN components using the empirical SED templates 


11 



Fig, 8. — UV-1R SED models of W2332— 5056. The dEished linra represent the different components in the SED, including starburst (grey), 
normal ditik (black)* and AGN (red and green), bi^ed on the empirical SED templates from Aasef et al. (2010). Because of the variability 
between the different epochs, (see §3 1), we separate the fitting into three different epochs; 1991 August - 1992 November (GALEA" -f- DSS 
- blue). 1999 October (GALEA A 2MASS - green), and 200S November - 2010 May (GALEX + + WISE - red). The DSS-B 

data (grey dot) taken in 1975 is not included in the fitting. The UV data from GALEX is fixed for fitting all epochs (see text for details). 
The earliest epoch data (blue solid line) is well fit by a galaxy without any AGN component. The later two epochs (green and red solid 
lines) are fit by adding an AGN component which brightens over time. 


of AsseF el uL (2010), C)nc challenge is that the source 
is clearly variable at optical wavelengths (see §3.1 and 
Figure 0), Indeed, the red optical SED observed by Mag- 
ellan in 2008 is difficult to reconcile with the very blue 
z* — J color implied by the 2MASS ./-band photome- 
try obtained in 1999 unless we invoke variability. We 
therefore divide the photometric data into three epochs 
{1991.7 1992.9, 1999*9, and 2008.9-2010*5, in roughly 3- 
year-long time bins) for SED lifting purposes (Figure s). 
The multi-epoch SED fitting was done requiring a com- 
mon, non- variable host galaxy using a combination of 
empirical star burst, disk galaxy, and elliptical galaxy 
Templates from Asset' et ah (2010) to represent the differ- 
ent age stellar populations in \Y 2332—5050* We required 
these SED galaxy components to match t he morphologi- 
cal modeling discussed in §3.5; namely, we required that 
the r* fluxes from the three galaxy components sum to 
one- third of the total r f flux, or ~ 60/j.Iy. For the AGN 
component, we adopted a modified AGN SED template 
in which the warm dust emission at A > 1/nn is replaced 
by a power-law emission to coutinue the accretion disk 
emission to longer wavelengths. This modification was 
made in order to satisfy a reasonable fit between r f band 
and W2. We also apply foreground dust, extinction to- 
ward the AGN to accommodate the red optical color. 
This modified obscured AGN component is reminiscent 
of the synchrotron emission t hat characterizes blazars at 
these wavelengths (Abdo et ah 2(11(1); the WISE colors uf 
W2332— 5050 put it. close to the blazar strip of Massaro 
et ah (2012), Blazars, conventionally considered to he 
the radio galaxies with relativistic jets close to the line 
to sight, exhibit large amplitude luminosity variations 
over hours to days (see review by Ulrich et ah 1997)* 


However, the interpretation of W2332— 5050 as a blazar 
is inconsistent with the radio morphology", nor do we see 
strong evidence for short-time scale variability at optical 
and inid-lR wavelengths* We do not address this blazar 
interpretation further here* 

The results of the model fitting are presented in Fig- 
ure s, The 6VI L EX obs i * r xi 1 1. io i is , ob t ai nc id be t weei i 2 0 03 
and 2005, would appear to present a challenge. On 
the one hand, they could be associated with a variable 
AGN component observed during an epoch when we do 
not have other photometry available! to use' in the SED 
modeling. On the other hand, they could be associat ed 
with the non- variable host galaxy, or a combination of 
the host and variable? AGN. In practice, our results do 
not depend strongly on the GALEX data. If we com- 
pletely ignore those UV data, our best, fit model implies 
a heavily obscured AGN during the latest, epoch with 
E(B — V) = 0*81 ±0*05* For this model, t he host galaxy, 
dominated by the starburst (e.g., irregular galaxy) com- 
ponent, goes through the GALEX photometric points 
with the best fit luminosity to the irregular component 
being T jin ~ {2*fi ± 1.5) x 1O 1O L 0 , Additionally, the 
xlcrf increases from 1.G5 (dof = 6) to 3*2t> (dof — 7) 
if we exclude the irregular component in tlie SED fit- 
ting. These suggest that the GALEX photometry cor- 
responds to the host galaxy with little AGN contribu- 
tion* If we redo the fitting under that assumption, the 
resulting AGN obscuration basically remains the same, 
E(B -V) = 0,82 ± 0*0 1, and the starburst. component lu- 
minosity becomes slight ly more robust, but is essentially 
unchanged* L\ m — (3.2 ± 1*0) x IQ 10 X, . It is this latter 
t ■ ase w he re vvi ■ ass ume the dec ade- old GA L EX photom- 
etry is due to the host galaxy* which we show in Figure 
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S. 

An obscured AGN dominates the optica] and WISE 
photometry for the most recent epoch, while ten years 
ago, when the 2MASS photometry was obtained, mtr 
modeling is based on fewer bands but suggests a slighi ly 
less luminous AGN with a similar level of obscuratimu 
E{B — V) — 0.8. The AGN reddening in these epochs 
is slightly higher than t ho E(B - V) 0,4 ±0.1 implied 
from the X-ray column, Nh ** 2,3±Q.5x 10 21 cm” 2 (§3. t), 
if the conventional Milky Way value of Njj/E(B - VO ~ 
5,9 x 10 21 cm” 2 (Boldin rf. al. 1978; Ruehford et ah 21)02) 
applies. Finally, the same fitting was applied to the DSS 
photometry for epoch 1991.7 1992.9. Remarkably, I he 
best-fit model show's no AGN emission component is nec- 
essary, and is relatively well modeled by the host galaxy 
component alone. The DSS 5-band, taken in 197f», is 
not included in the fitting due to its early observation 
epoch, but appears not to affect the fitting results if in- 
cluded. The best-fit AGN amplitude is nominally zero, 
suggesting that only ^ 20 years in the past W2332— 5956 
did not have appreciable nuclear activity. We note that, 
although the best tit model has relatively larger uncer- 
tainty due to the quality of the DSS photometry, the 
strong variation in the optical band over two decades is 
significant. The interpolation of our high-quality optical 
photometry measurements in the current epoch (2008 
2010) suggests a ilux of ^ 300 ^Jy in the DSS /-band 
measurement, > lGcr higher than the measurement of 
SO ± 20 /iJy in 1991-1992 epoch. 

3.7. Host G a last i a PRONGS Candidate 

FR-I1 radio sources are almost invariably found to hr* 
associated with large elliptical galaxies (Matthews et al, 
1904; Owen A; Laing 1989). Such hosts typically have de 
Vancouleiirs 1 r 1 ^ 4 profiles (Owen A: White 1991), corre- 
sponding to a Sersic index of 4. However, a class of radio 
sources hosted in star forming galaxies has been identi- 
fied by Norris et ah (201 Hi), and later named as Powerful 
Radio Objects Nested in Galaxies with Star formation 
(PRONGS; Man et. al. 2010; Norris et ah 2012), These 
objects are radio-loud AGNs {> 10 24 WHz" 1 ) with ra- 
dio lobes < 200 kpc that are hosted by spiral or starburst 
galaxies. Mao et id. (2010) discovered 50 PRONGS in a 
systematic search of a 7 deg 2 field in tile Aust ralia Tele- 
scope Large Area Survey (ATLAS; Norris et ah 2000; 
Middclberg et ah 2008). Prior to this work, two radio- 
loud edge-on disk galaxies (NGC612 and 0313-192) were 
known in local Universe (Ekers et id. 1978; Ledlow et ah 
2001; Keel et ah 2006; Emontsct ah 2008). NGC012 has 
a collimated FR-l-like jet, and an FR-lMike lobe with a 
bright hot spot on the other end, sometimes described 
as a hybrid-morphology radio source (HYMORS; Gopal- 
Krislma k Wiita 2000), while (1313-192 has a typical FIM 
morphology. There are a few other similar nearby ex- 
amples where radio- loud objects are hosted in disk-like 
galaxies such as 30305 (Heckman et al. 1982), 3C293 
(van Breugcl et ah 1984), B2 0722+3(1 (Eti units et ah 
2009), PKS 1814-636 ( Morgan ti et ah 2011), and SDSS 
J 1409-0302 (Hota et ah 2H 1 1). They could be associ- 
ated with PRONGS but perhaps at a slightly different 
evolutionary phase. Previous work has suggested that 
PRONGS may represent an early evolutionary stage for 
radio-loud galaxies. An example is given by Norris et 


ah (2912), who present a YLBI image of a highly em- 
bedded, compact core-jet AGN with 1,7 kpc-long jets in 
IRAS F00 183-7 111, an ULIRG with L bo , - 9 x 10 l2 £ G 
and substantial star formation activity. 

The optical morphology (§3.5) suggests that the host 
galaxy of W2332— 5056 has a modest disk (Sersic index ^ 
0.9), not much bigger t han the Milky Way in size and lu- 
minosity. With a 5GJI z-to-4400A-corUin u urn flux density 
ratio > 10, the W2 33 2— 5056 system is clearly hosting 
a radio-loud AGN at wording to the definition of Keller- 
nianti el ah (1989). Although hosting a HYMORS-like 
source, the galaxy of the W2332— 505G system seems to 
be a normal disk galaxy, making it a candidate for the 
rare PRONGS classification. 

Based on its rest-frame NUV luminosity vL v ^ 1.1 x 
l() u ‘ /,. from the host galaxy SED, we estimate a star for- 
mation rate of 5 A/, , . y r 1 in W2332— 5056, using Et illa- 
tion 1 of Kennicutt (1998) with no significant dust atten- 
uation assumed. The total luminosity of W 2332—5056 
up to 22 fi ni is 3.7 x 10 n L . The non-detection by 
AKARI f IRAS , and SPT limit the FIR luminosity (±fir) 
to 3.5 x 10 12 L-:, if an Arp220-like FIR SED is assumed. 
However, this Arp22()-like SED model does m>t match 
the observed optical- to-near- 1 R photometry of the sys- 
tem. If a fraction of the star formation activity is highly 
obscured at optical wavelengths (Ay > 49) but releasing 
energy via dust emission, the 22/un (W4) flux density 
limits Lfui to < 3,5 x 10 11 i 0? or < 30 A/^yr” 1 (e.g,, 
10% that of Arp220; Ammtlnirauiuinli et ah 2000). Thus, 
the host galaxy of W2332— 5056 is not likely to be un- 
dergoing a major star burst. This star formation rate is 
c lose to a typical quiescent disk galaxy, unlike the basts 
of ot her PRONGS. 

sccti on Discuss! on 

Observation ally, as discussed in §3, W2332— 5056 
shows (i) unusual, complex radio morphology, including 
curved jets, and an emission complex within the central 
25 kpc, (ii) a disk-like galaxy as the host, of an FR-l/FR- 
11 hybrid radio source, (ill) long-term (3-20 yrs) optical- 
near-IR nuclear variability, (iv) DPE profile In H/3 and 
Ha, The standard interpretation of DPEs iu which the 
BLs originate from the surface of a geometrically thin 
accretion disk around a single central SMBH is not ruled 
out by the observed optical spectrum, and do not pro- 
hibit the observed photometric variability. Our best-fit 
single AGN model (§3.2.2) entails an inclined elliptical 
(e “ 9,6) thin disk ut a viewing angle of ^ 42° from 
the face-on. The lifetime of liigh eccentricity could be 
a concern, though the circularization of the disk is from 
inside-out; thus the outer disk can retain its ellipticity 
for longer periods of time (Sycr k Clarke 1992). On the 
other hand, a warped disk, or a single-arm spiral disk 
can mimic the BL profile similar to that from elliptical 
thin disk. These possibilities can be examined through 
the BL profile variability from a long term spectroscopic 
monitoring, and detailed line profile modeling, although 
there is no detectable variation in the line profile over 
the ^ 6 month period separating the two CMOS obser- 
vations. 

The radio emission features, however, indicate that 
W2332 -5056 has a more complicated configuration than 
a typical disk-emitter. First, the viewing angles of the 
primary jets and the accretion disk in the thin disk model 
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are not consistent with each other. The flux ratio of the 
two primary radio jets along the inner 4?5-9i # 5 is 2. 
This fairly small ratio suggests the primary jets are close 
to be in the face-on plane (the angle between the jet flow 
direction and the line of sight > 80°) unless the two jets 
are intrinsically different. The inclination of the best- 
fitted elliptical disk, however, is 30° -50° after relaxing 
other fitting parameters such as eccentricity to e = 0.2 - 
0.8. Secondly, the appearance of the radio emission com- 
plex (such as component U D” and “F” in Figure 3) within 
projected 25 kpc from the core suggests some event or 
continuous activities within the past 70,000 years if the 
complex originated from the core and propagates rela- 
tivist ically. Since the complex has high radio luminosity 
(/v/ IQ 40 ergs -1 ), and it does not have optical or mid-IR 
counterparts, the assumption that it originates from the 
AGN is viable. The causes of the radio complex could 
easily create imprints (such as a hot spot or a gap) in 
the geometric configuration of the thin disk, if the disk 
configuration is not part of the cause itself. 

We should note that if a rotating BH Ls misaligned 
with the accretion disk, the material will be realigned 
due to Lense-Thirring drag (Bardeen & Petterson 1975). 
SMBH-accretion disk realignment, which is sometime 
attributed to the cause of the XRGs (see §1), has a 
timescale below Myr for a SMBH of 10 s A/© under typ- 
ical AGN accretion rates (Dennett-Thorpe ct ul, 2002). 
We crudely estimate the SMBH mass of W2332-505G 
(with radio luminosity of — 5.7 x 10 40 ergs -1 at 5.6 
GHz and a core X-ray luminosity of L { 2 — lOkeV) = 
4.1 x 1G 44 ergs -1 ) to be ^ 4 x 10 s A/© using the funda- 
mental plane of black hole accretion (Merloni lit al. 2003; 
Gultekiu et al, 2009). The uncertainty on this mass is 
a factor of ^ 8 due to the scatter of the empirical ns 
lations between X-ray luminosity, radio luminosity, and 
BH mass. The large uncertainty in the SMBH mass will 
not affect the realignment timescale too much since it has 
a low dependency (oc to the BH mass. It is vi- 

able to consider that the radio complex is the result of jet 
reorientation due to the BH-disk realignment in the past 
0.1 Myr, although the cause of the initial misalignment 
between SMBH spin and the accretion disk Ls unknown. 
Under this scenario, the double-peaked BL profile is pro- 
duced by the warped disk induced in the realignment 
process, and the two curved radio jets W2332— 5056 are 
the relic of the early realignment stage. 

The other plausible scenario involves two SMBHs re- 
siding in the disk-like host of W2332-5056 as a result 
of a major galaxy-galaxy merger. Based on numerical 
simulations, disk formation takes ^1-2 Gyr after close 
encounter of the major merging event of the parent gas- 
rich galaxy pair (Robertson efc al. 2006). This period of 
time is enough for two SMBHs to sink down to the cen- 
ter of the merged system (Bcgebmm et al, 1989), and 
reach the stalling separation at a distance of order pc 
or less. In this case, double-peaked Balrner lines arise 
from a long-lived eccentricity or single-arm spiral disk 
induced by the dynamical interaction with the second 
SMBH within a misaligned orbital plane, as discussed by 
Eradeous et al. (1995); Bogdanovic et al. (2008). The bi- 
nary interaction could be responsible for disk warping or 
a precession of the plane of the disk, as indicated by the 
mismatched inclination angles of the disk between now 


(implied by the spectroscopic models) and a few hundred 
thousand years ago (implied by the close-to-equal bright- 
ness of the primary jets on the two sides), and leaves the 
curved radio jets (i,e, jet precession; Roes 1988) as well 
as the enigmatic radio features. 

This scenario might apply not only in W2332-5056, 
but also in some DPEs. A recent VLBI survey of SDSS 
DPEs with FIRST detections shows that 2 out of 6 
studied DPEs have compact radio cores at kilo- parsec 
scale separations (R. Deane, in prep.). The radio galaxy 
4C37.11 contains an SMBH binary system with 7.3 pc 
separation and ~ 1.5 x 10 s A/© in mass (Rodriguez et al. 
2006), In addition, there arc several radio-bright DPEs 
in which double nuclei with separations of a few kpc 
are seen in high resolution radio mapping, such as SDSS 
J 1536+0441 (Wrobel A: Laor2009; Bondi Peroz-Torrcs 
2010), SDSS J1502+1115 (Fu et al. 2011b), and SDSS 
J 1425+3231 (Frey et al. 2012). These results indicate 
that some DPE systems indeed host dual AG Ns. 

The current resolution limits of radio and millime- 
ter interferometers, such as VLBI networks and ALMA, 
reach ^ 5 mas (Middelberg & Buell 2008; Wootten & 
Thompson 2009), At this resolution, compact radio/mm 
sources with projected distances > 25 pc at z — 0,3447 
can be resolved. Although this high contrast aperture 
synthesis imaging can recover liigh surface brightness 
structures such as highly collimated jets near their launch 
point, 25 pc physical resolution can not resolve an SMBH 
binary at “intermediate” merging stage with the stalling 
separations of 0,01-1 pc (Begclman et a!. 1980; Yu 2002), 
Current resolutions can, however, study an early stage 
SMBH merger, when the scattering of circnmnuclear 
stars is the dominant process for losing angular momen- 
tum, 

4. SUMMARY 

We report the discovery of W2332-505G, a radio 
galaxy at z ~ 0.3447 with extraordinary radio morphol- 
ogy and optica] spectroscopic features. We summarize 
the observational facts as follows: 

(1) Exhibiting hybrid features between FR-1 and FR-II 
morphology, the radio continuum map shows a winding 
Doppler-boosted FR-II jet, a pair of curved jets, a bright 
core, and complex emission within 25 kpc. 

(2) The host galaxy of W2332-5056, unlike the elliptical 
hosts of most FR-II objects, is disk-like. 

(3) We see long-term 3-20 yr optical-near-IR brightness 
variations from the AGN component. 

(4) The optical spectrum of W2332-505G shows offsets 
of ^ 3800 km s - 1 between the centroids of the broad and 
narrow line components of Ho and 

We consider the following two scenarios to explain 
these unusual properties: 

(1) W2332-5056 has a single spinning SMBH misaligned 
with its accretion disk. The realignment of the accre- 
tion disk due to dynamical drag produces the broad- 
band Balrner profiles, leaving reorientated jets in the pro- 
cess, and making the complicated radio morphology. The 
cause of the initial misalignment is unknown. 

(2) W2332-5056 has a SMBH binary with a separation 
on the order of a pc. The dynamical interaction between 
two SMBHs induces the warping or spiral pattern in the 
disk of one SMBH which generates the broad line profile. 
The curved jets and radio complex within the central 25 
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kpc reflect some perturbation from the second SMBH on 
the accretion disk. 

Our results and analysis can neither rule out nor con- 
firm either of these two scenarios. The real nature of 
W2332-5056 remains a mystery. Further investigations 
will be required to resolve the complicated structure 
within 5" radius of the core and confirm their nature, 
and even probe the sub-components in the core if there is 
any- We are currently pursuing the high-resolution radio 
mapping of W2332-5056 to investigate our hypotheses 
of the radio morphology formation, and also monitor its 
optical spectrum and photometry. If our hypothesis is 
confirmed, W2332— 5050 will be an exciting system to 
study the physics of an SMBH merger event. 
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